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Abstract 

We study theoretically the physical properties of a magnetic impurity in graphene. The theory 
is based on the Anderson model with a very strong Coulomb interaction on the impurity. We start 
from the Slave-Boson method and introduce a topological picture consisting of a degree of a map 
and a winding number (WN) to analyze the phase shift and the occupation on the impurity. The 
occupation is linked to the WN. For a generic normal metal we find a fractional WN. In contrast, 
the winding is accelerated by the relativistic dispersion of graphene at half-filling in which case an 
integer occupation is realized. We show that the renormalization that shifts the impurity level is 
insufficient to invert the sign of the energy level. Consequently, the state at half-filling is stable 
unless a gate voltage is tuned such that the Fermi level touches the edge of the broadened impurity 
level. Only in this case the zero field susceptibility is finite and shows a pronounced peak structure 
when scanning the gate voltage. 

PACS numbers: 75.20.Hr, 72.15.Qm, 71.55.-i, 81.05.ue, 85.75.-d 
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I. INTRODUCTION 



Graphene has attracted many theoretical and experimental researches due to its unique 
properties which are also of relevance for technological applications l|, |2J. The hallmark of 
graphene, a monolayer of carbon atoms, is its electronic structure with the valence and the 
conduction bands touching at two inequivalent points K± at the corners of the first Brillouin 
zone (FBZ). The low energy dispersion around K± is relativistic (linear in momentum), with 
a massless Dirac fermion behavior [2J . Of a particular interest is the issue of how the nature 



of graphene is manifested in the behavior of magnetic adatoms 



3Ml7|| , a topic at the heart of 



many-body physics, and especially the Kondo effect The Kondo model with a linear 

dispersion and the Anderson model in ci-wave superconductors li| have been investi- 
gated already. As detailed below however, for the particular case of magnetic adatoms on 
graphene some additional features emerge. 
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9j. This is only valid at temperatures T > T^, where TV is the Kondo temperature. 
The anisotropic single channel Kondo model [10] and the Anderson model [11] for infinite 
Coulomb correlation (U) were also considered. A Fermi liquid behaviour 11[ were con- 



cluded. In contrast, Ref. 



121 ] arrives at a two-channel Kondo in graphene due to the valley 



degeneracy of the Dirac electrons leading to an over-screening and thus to a non-Fermi- 
liquid-like ground state. In Ref. [15|, we conducted a detailed symmetry group analysis 
to clarify the appropriate physical model and highlighted the various relevant symmetries 
that are realized depending on whether (A) the adatom is above one carbon atom or (B) 
this atom is in the center of the honeycomb. The contributions from the two Dirac cones 
are mixed. For the case B we found generally a multi-channel, multi-flavor Kondo model. 
While for A we inferred a one-channel, two-flavor behavior. To identify the correct starting 
Hamiltonian a symmetry analysis is imperative. For example, the detailed symmetry anal- 
ysis in [l5[ for the A and B cases yields that the realized symmetry groups are C^ v and Cq v 
point groups. To be consistent with these symmetries, the eigenstates for pure graphene at 
Dirac points should be recombined as to reflect the modifications imparted by the impurity. 
As a consequence, a single half spin with zero orbital angular momentum is decoupled from 
graphene in case B (in contrast Ref. |2Q| ). 

In this work, we focus on the situation A and consider a two-flavor Anderson model with 
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a relativistic dispersion relation. Therefore, the effect of the gate voltage can be studied 



in a wide spectrum since the charge fluctuations are already taken into account [2l|. We 
introduce a geometrical picture in form of a winding number (WN) and a degree of a map, 
to analyze the occupation and the phase shift. It is shown that the occupation takes on the 
values or 1 when the bare level is respectively above or below the Fermi energy at half- 
filling. In the presence of a gate voltage, our analytical and the numerical calculations show 
that the state remains stable only when the edge of the broadened impurity level touches 
the Fermi energy where a nonzero susceptibility occurs. These findings are traced back 
mainly to the relativistic dispersion in graphene. In section II, an appropriate formulation 
is worked out, followed by the topological interpretation in the section III. In section IV, 
numerical illustrations are displayed for the half-filling case and beyond. In section V, the 
renormalization of the impurity level and the occupation dependence on the varying gate 
voltage are calculated. We conclude with a summary of this study. 



II. FRAMEWORK 



We start from the Anderson Hamiltonian 15] 



H — H g + i^hyb + H ira 

where the terms H g , H hyh , and H ira describe respectively graphene, the hybridization, and 
the impurity. The graphene Hamiltonian reads 15 1 



E 



where = hvpk, and vp is the Fermi velocity, s and o are valley and spin indices. k c is the 
cut-off momentum that sets the linear dispersion region. c s k a is an annihilation operator of 
the one electron state \ska). The Hamiltonian of the impurity is treated in the U — > oo 



limit, in a standard way 
conservation 

where 
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25( . We introduce a bosonic field b to guarantee the (Q) charge 
Q = tfb + n f = 1, 
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and f a is the annihilation operator of an electronic state on the impurity with spin a. Within 
the slave-boson (SB) model H reads then 

H = e n f + H g + H hyh + \(tfb - 1). 

Here the normalized impurity energy level is 

£o = £o + A. 

The hybridization Hamiltonian is formulated as 



t>o is the hybridization strength, and Qo is the area of a unit cell. As usual the bosonic 
field is assumed to be condensed at the ground state and is described by a renormalization 
number as 



A and ( are determined by minimizing the free energy which leads to the equations 




(1) 



(&t> = (b) = c. 



C 2 = 1 - (n/> 



and 





The Green's function associated with the impurity is 




where z = u + — Sq — C 2 Sq(w + ), co + = to + i5, 5 — > 



(3) 



The selfenergy is defined as 




and can be given analytically as 




(4) 



The density of state 



of graphene around the Dirac points reads 



p(u) 



f2 oj 



27T(hv F ) 2 ' 
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Defining the local density of states (LDOS) of the impurity as 



•A/" fa 



3 

2n' 



G 



we find 



(n f ) = 2 / f(u)M fa du, 



where f(u) is the Fermi function. Following 22], we derive 23] 

(n f ) 



1 f £F 
— 3 / d\nz. 

7T J-D 



Here we introduced 



and 



z = r exp i[ir/2 — Q(u), A, £ 2 )], 



9(w,A,C 2 ) = tan" 1 



a; - go - C 2 ^S C 



(5) 



(6) 



is the phase of z in (— 7r/2, 7r/2). Therefore, the occupation number reads 

(n f ) = ^- [ F d (26 (w, A, C 2 )) = deg(the map), 
which is the Friedel sum rule for an impurity on graphene. 



(7) 



(8) 



III. TOPOLOGICAL INTERPRETATIONS 



Now we introduce a new picture that allows a topological interpretation by extending the 
concept of the degree of a map and a WN of a closed curve to an open curve. In Eq. (jSj), 
20(cj, \,( 2 ) defines a map: O h- > P, where O and P are both ID manifolds. O stands for 
the ID energy region from — oo to D; and the manifold P is S 1 . The integral in (jSJ) can be 
viewed as a winding process by varying the source point p (stands for ep). Simultaneously, 
the image point p scans in the manifold P shown in Fig. (la). If p has a cyclic winding, 
usually an integer for the number of times that the manifold O covers the manifold P is 
produced and called a WN (similar to the case of a continuous map [26]). Note, our degree 
of the map 2G has the same topological meaning for non-integer WN, it indicates then that 
the winding process is not complete. We note further, is the connection in the image 
manifold P, and <90 / du is the curvature of this manifold. 
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FIG. 1: (color online) (a) Schematics of the map from ID manifold O to P. (b) The winding process 
for a normal metal. 

For a comparison, let us recall the same map for a normal metal with a constant DOS 
(u) E (—00, £f\), he. 



When p reaches Eq, the image point p attains G = 0, meaning a half winding of the manifold 
P. When p moves over Eq and approaches the Fermi level Ep, the image point p stops 
somewhere in the upper branch of the manifold P if 2B 7^ 7r (see the endpoint in Fig. 
(lb)). The winding is not completed so that the occupation number is not an integer but a 
fractional number. In the case that a mean-field SB method is applied and the real part of 
the selfenergy is ignored, the map reads 



If £ is not zero, Eq is renormalized to become Eq which is mapped onto = 0. The A in e 
gives rise to a shift of the impurity level, and ( 2 changes the winding velocity of the image 
point. 

For the graphene case, we define a function 



which is shown in Fig. (2a). The starting point of the winding is fixed at 20 = —n since 
T — > —00 as u> — > —00. In Fig. (2a), the solid horizontal line indicates the position of 
Eq so that the crossing points with J 7 indicate the zeros ( 20 = ) of the image points p 
and the corresponding pre-image points (source points), i.e. u*s of p. Counting the number 
of the pre-image points from the inverse map, generally, there are two points, i.e. u;*,^, 
residing closely to u — —D in two separated regions, I and II, with opposite countings of 



Q(u) = tan^ 1 ((a;-£o)/Ao). 




^,A,C 2 )=^-C 2 9f?£o(^ + ), 



(9) 
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FIG. 2: (color online) (a) Numerical graphs for J-(u), A, ( 2 ) and (b) a schematic of the winding 

2 

process for ( 2 = 1, £q > 0. In (c) Q 2 = 0, Sq < 0. The energy unit is D/ir. In (a) r/ = ^ [rw° F y 2 * s 
assumed to be 0.4. 

the degree of the map. The index of oj* (u^) is +1 (—1). The behavior of u = (Dirac 
point at half-filling) is particularly important for our analysis. As the imaginary part of the 
selfenergy diminishes, the end of the image point for u> = is determined by the relative 
positions of the function J-"(u; = 0, A, ( 2 ) and e$. When the former is larger (less) than the 
latter, the endpoint is tt (—it). The linear dispersion around the Dirac points accelerates 
the phase shift so that the winding process speeds up to the boundary of the manifold P. 

A schematic diagram of such a winding for ( 2 = 1 and Eq > is shown in Fig. (2b). While 
the first cycle is finished in region (I), in region (II) p comes back from +tt clockwise to "A" 
point, approximately crossing the zero of 20 (at u^)- In region (III), p starts approximately 
from "A" winding anticlockwise. However, it does not reach the zero of 2B and turns back 
to —7i to finish the winding. Hence, a zero winding is concluded for such a case, which is 
consistent with the condition of minimizing the free energy. Fig. (2c) shows the winding for 
zero ( 2 and Eq < where the WN is 1. We should note that zero ( does not mean that there 
is no effect from graphene to the impurity state. The charge fluctuation still renormalizes 
the impurity energy level via the parameter A in this case. 

It is readily shown that the occupation takes on only the values or 1 when eq > or 
< 0. For instance, starting initially from e < 0, A = and ( 2 = we obtain the occupation 
1. Thus, ( 2 = is stable against the renormalization process. Starting with Eq < 0, A = 
and ( 2 ytz 0, we find the occupation 1 after one step. This leads to a new ( 2 = 0. Therefore, 
we may conclude that the renormalization in graphene is quite different to that in normal 
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FIG. 3: (color online) Determination of A. The solid lines correspond to y\ = (eo — £o)/2?]-D — 
(1 + Vg/D). The dashed curves are calculated from In terms in Eq. (jlOp . The dotted lines are grid 
lines. The energy unit is D. The other parameters are eo = —0.3 D, r/ = 0.02. 

metals. When the bare level is below the Fermi level the renormalization effect is small. 



IV. NUMERICAL ILLUSTRATIONS AND THE CASE BEYOND HALF-FILLING 

Let us consider the effect of a finite gate voltage v g when the system is away from half- 
filling. To calculate the susceptibility, we introduce a homogenous, static external magnetic 
field and let it tend to zero at the end of calculations so that a zero-field susceptibility is 
obtained. The occupation number and the magnetization of the impurity read 

and 

where the phase CT is now spin-dependent containing the spin-dependent S 0(J , and 

£oa = e - vh, h = ^gfi B B, 

where g is the Lande factor, fiB is the Bohr magneton. To investigate the stability of ( 2 = 
at non-half-filling, it is crucial to determine the position of A. We find, for ( 2 = 0, 



S 



A 



7/ 



2p + ^ g ) + E^~oln^ 



2(D-v g )+j: (T e \n- 



(10) 



\£Oa+D\\v g ~£ 0a 

for fg + cr/i < 0, and > respectively. We solve A graphically, as shown in Fig. 3. The 
crossing points of the solid lines and the dashed lines deliver the solutions of A for a given 
gate voltage. The logarithmic terms show an interesting behaviour, the peaks being present 
at the positions of the gate voltages (i.e. the positions of the Fermi levels). This graph differs 
substantially from that discussed by Lacroix {^J for a normal metal. A good approximate 
solution of A is inferred by replacing the renormalized impurity level in Eq. ( FlOl) with Eq. 
As known, the renormalization to the energy level by the hybridization is unlikely to change 
the sign of the level. The occupation can not be changed by A alone. As a consequence, 
( 2 = is stable as long as io a < v g . For a finite ( 2 we resort to numerical calculations. We 
derive the susceptibility (higher orders in v g are ignored) 

vC 2 »b sgn(?j g(J ) (e„ - r]( 2 v ga ) 



X 



(ii) 



where v ga = v g + ah. Interestingly \ changes its sign in accordance with the sign-change of 
v ga which reflects the particle-hole nature of the graphene. 

Figs. (4a) -(4c) show self-consistent numerical calculations demonstrating our above ar- 
guments. In Fig. (4a) A changes only slightly for v g > —0.1, meaning the charge fluctuation 
is not large when the Fermi level does not reach the impurity level. A increases when a suffi- 
cient negative gate voltage is applied. However, its value can not convert the bare impurity 
level from negative to positive (in the local moment regime in our study). In Fig. (4b), 
the occupation number varies with v g . For a small negative gate voltage, the fully occupied 
state is still stable, a delta-function type DOS is induced (since ( 2 ~ 0) which is schemat- 
ically shown by the right insert in Fig. (4b). When the Fermi level touches the impurity 
level, the charge fluctuation has a strong influence leading to a remarkable decrease in the 
occupation. A broadening of the impurity level occurs (linear in v g ). This strong variation 

the difference to our case 



in the occupation also shows up in (i-wave superconductors p 
is that the steep decrease stems from the full occupation under insufficient gate voltage in 
graphene. For a comparison, in normal metals, the occupation is not complete even before 



the Fermi energy touches the impurity level, leading to a much smoother change 
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FIG. 4: (color online) (a) Numerically calculated A, (b) (re/), and (c) zero field susceptibility 
Xh-^o- (d) A graphic method to obtain the occupation by setting A = 0.135. The energy unit 
is D/tt, £q = —0.34, 7] = 0.02. The shaded regions of the insets in (b) indicate the broadened 
impurity level and the vertical lines indicate the positions of the Fermi level with a gate voltage. 
The left (right) inset corresponds to lower (higher) occupation. The function in (d) is y\ = 
sgn(w g ) - £o/>g| + (1 - (n f ))r]sgn(v g ) In |vr 2 - v^/v^. 

It is instructive to determine the special occupations for the varying narrow region. When 
= £ r 0o- + C 2 ^So-> ( n f) = 0-5. When the gate voltage is lowered further, (n/) = 0.5+ e< , 

sgn(v grT ) i \D 2 -Vg 



where Q^(v g ) = tan 1 (w(v g )) and 



w(v e 



^s^ln 1 v2 Kl . In Fig. (4b), at this point, 



(rif) ~ 0.4. The Xh~>o is only nonzero when full occupation, i.e. (nf) = 1, is violated which 
is shown in Fig. (4c) by a peak as lowering the gate voltage. To understand this strong 
change in the occupation, we derived the solutions of the occupation by the graphic method 
shown in Fig. (4d). The crossing points between the dashed and dash-dot lines delivers 
the solutions of the occupation. The tangent function is deformed by the relativestic linear 
dispersion to a step-like function resulting in a steep change with v g of the occupation, which 
is also comprehensible from our winding picture. Since ( 2 = fixes the endpoint to 20 = n, 
the occupation can change only when the Fermi level crosses the renormalized impurity level 
where a sign change occurs. 
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V. DISCUSSIONS AND INTERPRETATIONS 



For a deeper insight into the step-like variation of the occupation with the gate voltage, 
we consider the velocity of the parameters A and (rif) with respect to e F in absence of an 
external magnetic field. A is governed by the relation 

-^J- D w+-z -C*z - (12) 

After some algebra we find that 

d\ _ (go - e F )Mf{e F ) 

de F ~ e ' [ ' 

where A// = XL-^/o-- The variation of the occupation with respect to e F can be derived as 

|i>, % H 9e| f >( ' ) (14) 
oe F oe F 

From Eq. (JHJ), we infer that the velocity of the occupation with a varying gate voltage 



is determined by the LDOS at the Fermi level. This velocity or LDOS also describes the 
curvature of the manifold at the Fermi level, as interpreted in the previous section. By 
noting the fact that the LDOS is positive, the occupation increases with raising the gate 
voltage above the Dirac point. When lowering the gate voltage below the Dirac point, the 
occupation decreases. To know how fast the velocity of variation can be, we write explicitly 

MM = 4( W -^o-CW + (C 2 9So) 2 - (15) 
When ( 2 — > ((n/) — > 1), the LDOS develops a delta function at the virtual impurity level. 
This is the case of half-filling. Therefore, if the virtual level is below the Fermi energy and 
a gate voltage is applied to drive the system away from the half-filling regime, the vanishing 
LDOS is manifested as a vanishing velocity of the occupation with the gate voltage, unless 
the Fermi level touches the virtual impurity energy. This is the reason for the behaviour 
observed in the numerical calculations. The variation of A with the gate voltage, From Eqs. 
( fl3|) and ( fl5|) . shows a maximum at = e F and the restriction of ( 2 disappears. When the 
gate voltage is apart from the virtual level, A decreases. This can also be observed in Fig. 
(4a). 
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VI. SUMMARY 



In summary, we investigated the Anderson model for a magnetic adatom above one carbon 
atom of one monolayer of graphene. We utilized a topological method, i.e. a degree of a 
map and a winding number, to analyze the occupation of the impurity and the phase shift. 
It is found that the phase shift is accelerated by the relativistic dispersion of graphene to 
complete one winding or zero winding when the impurity level is respectively below or above 
the Fermi energy at half-filling. The occupation varies dramatically from 1 (full occupation 
of the impurity) in a narrow range giving rise to a peak in the zero field susceptibility. The 
velocities of the renormalization of the impurity level and the occupation with respect to the 
gate voltage are worked out and an interpretation of the step-like variation of the occupation 
is provided and is consistent with the topological picture. 

The work is supported by DFG and the state Saxony-Anhalt. 
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